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ABSTRACT 


In an angle-measuring navigation satellite system using a 
pair of crossed interferometers located on a satellite in synchronous 
orbit, three parameters are needed to determine a user's position 
unambiguously: the phase difference between received signals which 

had been transmitted by the two antennas on each of the two inter- 
ferometers, and the user's altitude. 

The two phase difference measurements yield a line of possible 
user locations, and the addition of the altitude measurement reduces 
this line to a single point. Instead of measuring altitude via a 
barometric altimeter, a method is proposed here which makes use of 
the navigation signals received after reflection off the earth's 
surface. The iterative procedure used here employs the arrival time 
difference between direct and reflected signals. 

Based on previous calculations of errors in measuring the 
electrical parameters it is concluded that, for North Atlantic 
coverage and specular reflection, altitude measurements can be made 
with a l-o error of 65 meters. 
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1.0 INTRODUCTION 


In an angle -measuring navigation satellite system using a 
pair of crossed interferometers located on a satellite in synchronous 
orbit, three parameters are needed to determine a user's location 
unambiguously: the phase difference between received signals which 

had been transmitted by the two antennas on each of the two inter- 
ferometers, and the user's altitude. The satellite is assumed in this 
study to be transmitting at L-band and to be nominally located over 
the equator at 30°W longitude, thereby providing complete North 
Atlantic coverage. 

Each of the two phase difference measurements places the user 
on the surface of a cone whose vertex is the satellite (an approxima- 
tion good for great distances from the satellite) . The intersection 
of these two cones is a straight line extending from the satellite 
towards the user. 

A knowledge of the user's height above the earth reduces this 
line of possible positions to a single point. Normally, the altitude 
would be determined by a barometric altimeter. An alternate altitude 
determination method, which is proposed here, makes use of the satellite 
signals received after reflection off the earth's surface. This pro- 
cedure, which involves measuring the arrival time difference between 
direct and reflected signals, could be used as a backup for the 
barometric altimeter, or even in place of it. 

Several types of signaling (AM, PM, RF pulses) have been 
discussed in the literature (Ref. 1, 2, 3) and will not be repeated 
here. It will be mentioned, however, that regardless of the signaling 
employed, separation of the direct and reflected signals is necessary 
for accurate measurement of the arrival time difference. 

In this study it is merely assumed that the user is able to 
make the proper phase difference measurements and time delay measure- 
ments, regardless of the type of signals used in the system. 

A brief description of the geometry of position-finding in an 
angle -measuring system is found in Sec. 2.0. Section 3-0 is concerned 
with an analysis of the height determining procedure which makes use 
of multipath information, while Section 4.0 deals with an analysis of 
the errors involved. The results of a system evaluation and some 
conclusions are found in Section 5-0. 
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2 . 0 MATHEMATICS OF ANGLE MEASURING NAVIGATION SATELLITE SYSTEM 

The synchronous satellite will be considered located nominally 
over the equator at 30°W longitude. One interferometer boom is located 
in the equatorial plane in an east-west position, while the other inter- 
ferometer boom is in a north-south orientation. A nominal coordinate 
system, x 0 - y 0 - z Q , is considered with the origin at the nominal 
satellite location. This nominal coordinate system is shewn in Fig. 1. 

The satellite will usually be at some location other than its 
nominal one, and for this reason it is given a local coordinate system, 
x-y-z, also shown in Fig. 1. 

The satellite altitude is constantly monitored by ground-based 
tracking stations which relay the pertinent information to the field 
of users at regular intervals. 

The matrices used in transforming the nominal coordinate system 
into the local coordinate system, and vice versa, are found in the 
Appendix. 

Consider the pair of antennas along the x-axis, as shown in 
Fig. 2. The measured phase difference at the aircraft between the 
signals transmitted by each antenna is given by 

Acp = 2n ^ cos tlr - n2« (2-1) 

T x A. *x 

where A. is the wavelength of the L-band signal (1.6 GHz), and n is 
the largest integer in D/A. cos *jr x , since measurements are made modulo 
2n. Clearly, there are several values of \|r x which would yield a 
particular measured value of A<P : 

♦x = oos ’ 1 [jTd n2,,) ] • (2 ' 2) 

A user, then, would be located on the surface of a cone whose 
vertex is the x-y-z origin, whose axis is the x-axis , and whose half- 
angle is i|r x . Each cone, corresponding to each possible value of i|r , 
intersects the earth forming a curve. It can be shown (ref. l) that 
if the boom length, D, is 20 wavelengths, the minimum distance between 
these curves is 1000 nautical miles. Thus, if a user knows his posi- 
tion to within 1000 nmi, he then knows the proper cone and proper value 

of A . 

T x 


The equation of the cone is given by 
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Fig- 1 General Satellite Coordinate System 
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Fig. 2 x-Axis Interferometer Geometry 
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(2-3) 


,12 2 2 
♦ x x - y - * = °- 

Actually, eq. (2-3) describes a pair of cones, the second one of which 
is shown in Fig. 2 by the dotted lines. The presence of two cones is 
of importance only when their intersections with the earth are close. 
However, since a user will be making new measurements every minute or 
so, for instance, he will be able to determine the proper choice of 
cones . 



Similarly, by measuring the phase difference between the 
signals from each of the two antennas on the y-axis, the user finds 
himself situated on the surface of a second cone with vertex at the 
x-y-z origin, with axis along the y-axis, and with a half-angle of 
\J/y This cone’s equation — again, there is really a pair of cones — 
is given by 

i|r I y- - x - z =0. (2-4) 



The intersection of the cones is a straight line which passes 
through both the satellite center and the aircraft, and which is 
described by the following equations: 


where 


x _ y = z 
a b c 


a = + cos i|r 
Y x 


b = cos \lr 

■ y 


(2-5) 


> (2-6) 


c = 


+ *\j 1 - COS 


cos 


J 


The reason for the choice of signs in eq. (2-6) is easily 
seen in Fig. 3> Since there are two cones centered on the x-axis and 
two centered on the y-axis, their intersections result in four straight 
lines. The field of users will be sent data regarding satellite atti- 
tude as determined by the ground stations, so that determination of 
the proper signs for a and c is easily obtained. 

The center of the earth (assumed to be a sphere of radius R e ) 
in the x-y-z system is at the point (x e , y e , z e ), and the equations 
of the line from the satellite to earth center are given by 
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(2-7) 



The distance, d g , from the satellite to earth center is given by 



(2-8) 


The line from the satellite to the aircraft intersects the 
earth at the point (x^, y^, z^), where 


X. = 

aT 


1 

*1 = 

bT 


z . = 

cT 


1 

T = 

s 

OJ 

CO 

S = 

ax + 
e 

by, 



+ cz 

e 




-\ 

> (2-9) 


J 


The derivation of eq. (2-9 ) is found in Appendix A. 

The angle, 6, between the satellite-earth center line and the 
satellite-aircraft line is given by 


9 


-1 

cos 


_S 

d 

e 


( 2 - 10 ) 


V 


Once the user's altitude is determined (see Sec. 3-0) and the 
point on earth directly beneath the aircraft, / x , y , z "N is 

V. °AE °AE °AE/ 

found in the nominal coordinate system, the aircraft's latitude and 
longitude are then computed by the following equations: 


Lat . = 


Long. = 


The point (x , 
°AE 


. -1 
sin 



( 2 - 11 ) 


He + z 

_ n o -.-1 S °AE 
30 + cos 


R g cos (Lat.) * X o AE <: °' 


(2-1 2) 


y , z ) is referred to as the sub-aircraft position. 
°AE °AE 
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3-0 HEIGHT DETERMINATION METHOD USING MULTIPATH 

The height determination method employed here is an iterative 
one which makes use of signals received by an aircraft after specular 
reflection off the earth's surface. Besides measuring the two phase 
differences discussed in Sec. 2.0, the user will also measure the time 
delay. At, between reception of direct signals and reflected signals. 

It will be assumed that the specular reflection point, M, is 
in the satellite-aircraft-earth center plane, P. This is not an 
unreasonable assumption, since the boom length, D, is only 20\ 

(3-75 meters at L-band) . 

It will be necessary to consider another satellite-centered 
coordinate system, the u-v system, shown in Fig. 4, such that the u-v 
plane is the same as the plane, P, and such that the u-axis is directed 
towards the earth's center and the v-axis is directed towards the 
v-coordinate of the aircraft. 

The details regarding the transformation from the x-y-z system 
to the u-v system, and vice versa, will be found in Appendix B. 

The sub-satellite point (u^g, v^) , or more precisely the 
point (x , y , z ), determines the aircraft's latitude and 

°AE °AE °AE 

longitude (eqs. 2-11 and 2-12). 

The height determining procedure is based on the following 
manipulations : 

Having determined values for 9 and At via the phase difference 
and time delay measurements, the user then determines the point (u^, v^), 
which is the intersection of the satellite-aircraft line with the 
earth, and imagines a tangent to the earth at this point, as shown in 
Fig. 5- Since this is the first of several points of tangency*, it will 
be denoted (u. , v. )• Using this tangent line, a satellite image point 
is found.** 1 1 1 1 


* It may be inefficient to start the iterative procedure with a plane 
tangent to (uj, , v.j, ) . However, the simulation of the system shows 
that only threi or four iterations are needed using this point; 
thus, it is not absolutely necessary to begin the procedure, with 
another point. 

** It is possible that the image method which is applicable to 
spheres would be useful. 
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The path length difference, A&> between the direct and reflected 
signal paths is related to the arrival time delay between them by the 
following equation: 



(3-1) 


where c is the speed of light. 

Consider , now , t he hy perbola of Fig. 6. The difference between 
the focal radii, P'F^ - P'F p , of any point, P' , on a hyperbola is a con- 
stant, and is equal to the length of the transverse axis, 1. One can 
therefore imagine a hyperbola with transverse axis length equal to the 
measured path length difference, A 4, as given by eq. (3-1). 

Consider, then, such a hyperbola, shown in Fig. 7, with foci 
at the satellite and satellite image point and with transverse axis 
equal to A The intersection of this hyperbola with the satellite- 

aircraft line is the point (u. , v. ) and is the first estimate of the 

_ 1 A 1 

aircraft's location. This point would be the true location, since it 
is on the satellite-aircraft line and since the observed path length 
difference at this point is Al, were it not for the fact that reflection 
takes place off the tangent plane (line) in this approximation and not 
off the earth's surface. 

A line is then imagined from the satellite through this planar 
point. This line intersects the earth at the point (uq p , vq ? ), as 
shown in Fig. J; a new line, tangent to the earth at this point, is 
considered, and the process is repeated. 

When the process has been taken through as many iterations as 
is necessary to satisfy required accuracies, the subaircraft point, 

( u pj?> V A]P 5 and aircraft's altitude, h, are then determined (Fig. k) 
by uhe following equations: 


h. = 
3 


^( d e ' u a3 2 + V 

^ n -I 


- R 


AE. 

3 


R u. + d h. 
e A e j 

si 

R + h. 
e 3 


AE. 

J 


R v. 
e A. 

si 

R + h. * 
e 3 


(3-2) 

(3-3) 

(3-4) 
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where the j refers to the values of the indicated parameters in the j th 
iteration. The details of this procedure are found in Appendix A. 

The point (u^g, v ae) ^ en converted back to the nominal 

coordinate system to become the point (xq^, Yopg’ z °ae^ which is 
used to determine the aircraft's latitude and longitude from equations 
(2-11) and (2-1 2). 
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'•1.0 ERROR ANALYSIS 


The altitude, latitude , and longitude of an aircraft may be 

written as functions of it , ilr , At, a, 8, v, a , b , and c . These 

. Y x Y y o o o 

measured parameters are independent, zero mean random variables with 

. 2 ' 2 2222 2 2 ,2 
variances a , a , a > a a , a , a , , c b and a c , 

T x T y o-o o 

respectively. 

The variance of altitude errors is then given by 


alt 


3 alt 

2 a 2 + 

3 alt“ 

2 2 

/"T 1 

3 

alt 

* V ♦ 

~3 alt’ 

3 \(f 
_ x 

Y x 

3 lit 

L y 

•a + 

Y y 

3 

At 

3 a 


a 


3 

alt 

2 2 

"3 

alt 

2 2 
•°Y + 

3 alt’ 

2 a 2 

3 alt 

3 

P 

‘ °3 + 

3 

Y 

3 a 

L o _ 

• a + 

a 

o 

3 b 

L o J 


“|2 


3 alt 


. a c 

L O. 


(4-1) 


and is a function of aircraft location. 

Similar expressions may be written for the variance of latitude 
errors and of longitude errors. 

The derivatives involved were all expressed analytically and 
were included in the computer program as part of the system evaluation. 
Eor instance, the altitude at the jth iteration is given by 


hj ” i (< V U A d >2+ V A. 2 - R e 


(4-2) 


and the partial derivative of h. with respect to, say a, is 


3h . 


3cx R + h . 

e a 


/ 3d g 3 u A.\ 

^ d e~ u Aj^ a " 3 a ) + V A.. 3 a 


9v a.~- 
1 


(4-3) 


Clearly, each of the three partial derivatives in eq. (4-3) now have to 
be evaluated. From 
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2 


(4-4) 


d 

e 




+ 


z 

e 


the partial derivative is 




9x 

ri + y e 



z 

e 



(4-5) 


and the partials in (4-5) now have to be evaluated, etc. The nine 
derivatives in (4-l),as well as those involved in the variance of both 
latitude and longitude errors, were all evaluated in this manner and the 
corresponding expressions made a part of the program which evaluated 
altitude, latitude and longitude. 


The following values were used in the system evaluation: 

a = a =50 prad [Ref. l] 

*x *y 

a , = 0.1 ns [Ref. 2] 

At 


°a = °3 = °y = 1,0 |irad ^ Ref " 


a a = cr^ = a c = 20 meters [Ref. 3] 
o o o 


It should be pointed out that not all possible error sources 
were considered in this study. The earth, for instance, was taken to 
be a sphere, and reflection was assumed specular. These factors should 
be considered in any further work on this subject. 
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5-0 RESULTS AMD CONCLUSIONS 

The system was evaluated by considering various values of the 
input parameters (i|r , At, a, 0, y, a Q , b , c ) corresponding to four 

^ y 

altitudes at each of the North Atlantic locations denoted in Fig. 8 by 
the letters A through H. These particular points were examined since 
the coverage area under consideration extends from 40°N to 6o°N and from 
10°W to 50°W. 

For each set of input quantities, the corresponding altitude, 
latitude, and longitude were computed, as well as the standard deviation 
of their errors. The results of these computations are given in 
Table 4-1. 

Figures 9 and 10 are based on the values given in Table 4-1. 

Figure 9 shows a. „ plotted against relative longitude (that is, longi- 
Alii 

tude relative to that of the nominal satellite position), while Fig. 10 

has a T . m and a T _™ vs relative longitude. 

LA. I LUNUr 

An examination of Table 4-1 and Figures 9 an d 10 yields the 
following conclusions regarding aircraft position determination using 
multipath: 

a. Altitude determination errors increase with increasing 
aircraft altitude and with increasing latitude. 

b. The greatest lcr altitude error was on the order of 65 meters 
and occurs at 6o°N at an altitude of 19 km. This figure is 
within the limits of a barometric altimeter (Ref. 2, 4). 

c. Latitude errors and longitude errors increase with increasing 
latitude and with increasing relative longitude. 

d. Position fixes can be made with la accuracies of 1*5-3 nmi 
in latitude and 1-1. 5 nmi in longitude in the North Atlantic 
coverage area. 

A few words should be said regarding the convergence of the 
iterative procedure. A method is described in Appendix C in which 
At, h, u A , and v. are determined geometrically with 9 and a known 

(see Fig. 4). This computed value of At is used in the iterative pro- 
cedure with the same 6. The values of h, u^g, and vae computed here 
agree (to at least l4 places) with those found geometrically for each 
North Atlantic position under consideration (see Fig. 8 and Table 4-1). 
For all extent and purposes, then, the procedure could be said to 
converge in the area of interest. 
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Fig. 8 Area of Interest in North Atlantic Showing Locations 
Considered in System Evaluation 


Table 4-1 


1-Sigma Errors 


Altitude 

(km) 

°ALT 

(m) 

1 

a LAT 

(nmi) 

a L0NG 

(nmi) 


18.967 

65 




14.440 

54 



A 

10.686 

49 

2.83 

1.50 


7-734 

47 




18.903 

6l 




14.484 

52 



B 

10.738 

47 

2.82 

1-31 


7-473 

45 




18.653 

56 



C 

14-377 

50 

2.78 

1.18 


10.840 

46 




7.665 

43 




18.852 

55 



D 

14.233 

49 

2.78 

1.09 


10.729 

45 




7.508 

4i 




18.929 

26 




14.222 

25 



E 

10.681 

25 

1-51 

1.18 


7-130 

24 




18.912 

25 




14-377 

24 

1.49 

1.11 

V 

IO.966 

24 




7-548 

23 

1 




18.856 

24 



G 

14.439 

23 

1.48 

1.06 


IO.674 

23 




7-346 

23 




18.918 

24 



H 

14.153 

23 

1.47 

1.03 


10.897 

22 




7-635 

22 
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(Meters) 



RELATIVE LONGITUDE (degrees) 


Fig. 9 1-Sigma Altitude Errors vs Relative Longitude 


LAT 1 59 . 7° 


= 40.4° 


20 




6.0 


REFERENCES 


1. Perry I. Klein, "An Analysis of Advanced Navigation Satellite 

Concepts," Fh.D. Dissertation, Univ. of Pa., May 1968. 

2. Pierre A. Portmann, "A Study of the Capabilities and Benefits 

of Navigation and Traffic Control Satellites - Volume III: 
Study of Altitude Determination from Satellite Signals," 
Stanford Research Institute, Final Report on NASA Contract 
NSR 05-019-243, Sept. 1970. 

3. RCA, "Phase Difference Navigation Satellite Study," Inter im 

Scientific Report on NASA Contract NAS-12 -509, June 1967. 

4. G. R. Thomas, "Analysis of a Spinning Interferometer Navigation 

Concept Study," Ph. D. Dissertation, Univ. of Pa., May 1971. 


22 


Appendix A 


DERIVATION OF EQUATIONS 


1. Intersection of 2 Cones (Satellite -Air craft Line) 


2 2 /, 2 \ 2 
y + z (tan ^ ) x =0 


2^2 2 , v 2 . 

x + z - (tan 4) y = 0 

O' 


(A-l) 

(A-2) 


Solving (A-l) and (A-2) for z and equating , 


2 2, x 2 2 (, 2 . v 2 2 

z = (tan i ) x - y = (tan i|f )y - x 

X j 


x 2 (l + tan 2 ) = y 2 (l + tan 2 \|r ) 
x y 


2 , 

COS Tjf 


2 , 
cos \|r 


- cos ilr 
T x 

From (A-2), x 2 = (tan 2 i|r )y 2 - z 2 

y 

Then (A-l) becomes 


COS l|( 


y 2 + z 2 - (tan 2 f x ) jjtan 2 i|r y )y 2 - z 2 J = 0 


y 2 + z 2 - (tan 2 i|r )(tan 2 ijf ) y 2 + (tan 2 i|f )z 2 = 0 

x y x 


(A-3) 


y 2 (l - tan 2 i|; • tan 2 \(f ) + z 2 (l + tan 2 i|i x ) = 0 
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Let 


(Note: 


y 


1 + tan ik 

]_x 

2 2 
1 - tan \1( • tan i i 
Y x Y y 


2 , 

cos \fr 


x 1 - 


sin^ilj sin^ilf 
Y x 

2 , 2 , 

COS f COS \Jf 


2 . 
cos \Jr 


2 2 2 2 
cos i|r cos * - sin t|r sin ik 

Y x Y y Y x Y y 


2 , 

cos \j; 


.2 .2 2 2 
sin sin ijf - cos \k cos ik 
Y x Y y Y x Y y 


2 2 2 2 
(1 - cos ^ x )(l - COS ljfy) - COS (jf x cos^iji^ 


2 2 2 
cos tk 1 - cos ik - cos \k 

y Y x Y y 


COS l(f 


= + 


r-. 2 ' 2 ~ 

\ 1 - cos ik - cos ik 
Y x Y j 


(A-4) 


cos ijf cos ij^ 


= + 


2 2 
Ml- cos \|/ x ~ cos i|r 


y 


(A-5) 


a = i cos i|f 


b = cos t(i 




c ~ 1 - cos2 ij/^- cos^ i|; 

a2 + t»2 + c^ = l) 


J 


(A-6) 
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Then 


x _ jr _ z 

a b c 


(A-7) 


is the set of equations for the satellite-aircraft line. 

2 • Satellite-Earth Center Line 

The equation of the line between (x, , y, , z ) and (x , y p , z ) 
is given by x l. x d d d 


x - x. 


x r x 2 


y - y-L 

hr y 2 


Z - z. 


V Z 2 


(A-8) 


Here, satellite: (0, 0, 0) 

earth center : (x . y , z ) 
e *e e ' 


x - 0 _ y - 0 


0 - x 


0 - v 


z - 0 
0 - z 


or 


x_ 

x 


= 2L. = 2- 

y, z 


(A-9) 


e " e e 

This, then, is the set of equations for the satellite-earth center line. 

3- Satellite-Earth Center Distance 

The equation for the distance between (x^ y^ z 1 ) and ( x y 2 , z^) 


is 


d e = 1 ( X x - ^ ) 2 + (y x - y 2 ) 2 + (z^ z 2 )‘ 


(A-10) 


d = M x + y ~ + z 
e \ e /e e 

4. Intersection of Satellite-Aircraft Line with Earth 

Equation of earth: 

(x - x e ) 2 + (y - y g ) 2 + (z - z g ) 2 = H g 2 


(A- 11) 


(A- 12) 
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Sat. -AC line: '- = £ = - (A-7) 

a b c 

Substituting (A-7) into (A-12) 

(x - x e ) 2 + ( I X - y e ) + (f X 

or (a 2 + b 2 + c 2 ) x 2 - | (a Xg + by g + cz g )x + (d g 2 - R g 2 ) =0 

a<_ 


- z ) =R 
e e 


x 2 - — (ax + by + cz )x + (d 2 - R 2 ) = 0 
2 a e e e e e 

o 


Solving and denoting the root by x^. 


x.=a.J (ax e+ by e+ cz e ) 


\ (ax g + by g + cz g ) 2 - (d g 2 - R g 2 ) 


Similarly., 

y i = b - 1 (ax e + b V oz e ) ' 


'l (a* e + by e + cz e ) 2 - (d/- R e 2 ) 


(A-13) 


t = <=■■{ ( ax e + t'y,,+ cz e ) - '! (ax 0 + by e + cz e ) 2 - (d e 2 - R e ) 


Let 


S = ax + by + cz 
e J e e 


T = S - 4 S 2 - (d g 2 - R g 2 ) 


(A-14) 


The point (x^, y^, z^) is the intersection (on the satellite side of 
the earth) of the satellite-aircraft line with the earth. 
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5 . Angle between Satellite-AC Line and Satellite-Earth Center Line 

X "V* z 

The angle, 9, between the line, — = f— = — , and the line, 

a l b l C 1 

x y z . , 

T~ = y- = — , given by 

a 2 b 2 C 2 


8 = cos 


-1 


a l a 2 + b l b 2 + C 1 C 2 


{ 


(A-15) 


2, . 2. 

'I + b l 


a,-+ b,“+ c x ^ & 2~ + b ^ + c ' 


2 . ’ 2 . 2 

i_ + c„ 

2 2 


8 = cos 


-1 


ax + by + cz 
e 'e e 


^1 a 2 + b 2 + c 2 x g 2 + y e 2 + z g 2 


(A-l6) 


or, from (A-6), (A-ll) and (A-lU) , 


„ -IS 

8 = cos — 
d 

e 


Latitude, Longitude 
From Fig. A-l 


sin a 


AE 


lat R 


(A-17) 


a, , = sin 
lat 


Also 


cos 


(? a long) 


,-xSe 

R 


\ _ cos a 

cos C a iat^ 


and 


cos a = 


R S + Z °AE 


R 


(A-18) 

(A-19) 

(A-20) 
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(The remainder of the derivations will use the general subscript, j • ) 


(b) Line tangent to earth at ^u^ , j 

Vj y 

S. = slope of curve (earth) at/u. , v^ j 

d - u. 
e 1 . 

_ 1 


v. 

l . 

3 


(A-26) 


as seen in Fig. A-3- 

The equation for the line tangent to the earth at that point 
is given by 


or 


v 


v 




v . 


(A-27) 


(A-28) 


(c) Coordinates of j satellite -image point 

From Fig. A-4 it is seen that the slope of the satellite- j 
satellite image line is equal to the slope of the j^h tangent point- 
earth center line and has the value 


v. 



Therefore the equation of the satellite-j satellite image line is 

(A-29) 



The intersection of the satellite-j satellite image line and the j 

tangent line occurs at (u c ., v c .). To find this point the equations 

J J 
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• u 


for the two lines (eqs. A-28 and A-29) are equated 

2 . 2 


d - u. 
e 1 . 

JL 

v. 

l . 

J 


/ 2 2 \ 
U. - ( U. + V. \ 
2 1 . ( 1 . 1 . ) 
j V j o / 


V . 

1 , 


u - 


V. 

1 . 
J 


d - u. 
e i , 


Solving for u. 


u c. = (d e - 
J 


V ' Z i 1 


and 


c . 
J 


v. • I. , 

x i 3 


(A-30) 


where 


/ 2^ 2 

d u. - ( u. + v. 

e l . V i . i . 

r = a ^ a_ 

3 (d - Uq .) 2 + V. 2 

v e 1 j / l. 


(A-31) 


line. 


Since the tangent line bisects the satellite-satellite image 


and 


u T - 2u 

I. c . 

J J 


1 


v T = 2v 
I . c . 

J J 


(A-32) 


or 


and 


U I. = 2 < d e" u U )T < 


y J 


2v. I. 


(A-33) 


(d) Coordinates of j estimate of aircraft location 


The angle, Yj> between the satellite-earth center line and the 
satellite-j^h satellite image line is given by 


v. 

Y . = tan" *4- 

0 d - u n - . 

e x j 


(A-34) 


as shown in Fig. A-4. 
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Consider a new coordinate system, u'- v', as shown in 
Fig. A-5- The u'-axis is along the satellite-j^h satellite image 
line. The transformations from the u-v system to the u'-v' system, 
and vice versa, are given by: 


cos y. - sin Y. \ , u 

J J ' 


sin Y . cos Y . ' x v 
J J 


(A-35) 


cos Y . sin Y . i 

J J \ / u' 


- sin y • cos Y • ' \ v* 

J «J 


At is the time delay between a signal received directly and a signal 
received by the aircraft after reflection. This is equivalent to a 
path length difference, frJi, where 



(A-36) 


Consider a hyberbola whose principal axis (the axis which 
joints the two foci) is the u'-axis, whose foci are the satellite 
and the j th satellite image, and whose center is (u' C j , 0) . Now, 
the difference between the focal radii of any point on a hyperbola 
is a constant, and is equal to the length of the transverse axis 
(the portion of the principal axis included between the vertices). 
Consequently, this length will be set equal to the required path 
length difference, above. 

The equation for this hyperbola is given by 




where 2Q 1 = fri, 





= 1 


(A-37) 


(A-38) 
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(A-35) 


and u' = u cos y. - v sin y. 

CJ C. J C. J 


• 0 = M 

• • ^1 2 


and 




(A-39) 


Thus, the equation for the hyperbola becomes 


u ■ - u 


■0* 


(vr 


= i. 




2 -(IT 


(A-40) 


The satellite-aircraft line is given in the u'- v' system by 


v' = tan(0 + y ) ■ u' 


(A-41) 


The intersection of the satellite-AC line with the 
hyperbola is the estimate of aircraft location, shown in Fig. A-5: 


u 


A. 

3 


or 


A. 

3 


A. 

3 


-(| At) 


u' c cos (0 + Yj) + 2 At 


cos(0 + y.) 

3 


v' = u' • tan(0 + y,); 


A. 

3 


c A , \2 
2 A ^ 


"A 


u' cos(e + y.) + t At 

c. 3 2 

J 


cos 8 


v A = U A ’ ^ an ® 

J J 


(A-42) 


(a-43) 
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th . , , .th 

The line between the j satellite image point and the j 

aircraft location estimate is given by 


u - u 


A. 

_i - 


V - V 


U A.“ U I. 


A. 


V V I. 

3 3 


or 


v = 


V A." V I. 

3 i u _ _3 — a 


U I. V A.“ U A. V I. 


3 ._I 


U A.“ U I. 
3 3 


U A. “ U ] 
3 


(A-44) 


Let 


m. 

3 


V A . ~ V I. 

__j a_ 

U A." U I. 

3 3 


and 


b . 
3 


W U A V I. 
3 3 3 3 

U A . ~ U I. 

3 3 


(a-45) 


Then, 


v = m .u - b . 
3 3 


(a-46) 


The intersection of this line with the j th earth tangent line (eq. A-28) 


is the point (u t ^. , v t . ) : 


u 


t . 
3 


/ 2 2 \ 
d u. - b .v. - ( u. + v. ) 

e V J ".1 V X 3 X 3 ) 


d - u. - m . v. 

6 X 3 J X 3 


v. - m.- u, - b . 

*3 J *3 ° 




(A-47) 


Let 


2 2 
W . = u , + v, 

■> h 4 


(A-48) 


The line from the satellite through (u^_ , ) is given by 

3 3 

V t. 


v = — “■ u 

■ u t. 

3 


(a-49) 
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The intersection of this line with the earth (see Fig. A-6) is the 
point /u. , v. \ : 

V Vi x j+i ' 


u. 

\j+l 


v. 

^+1 


u - aI R Sf.- d 2 V, 2 

lh I e 3 e h 

w. 

0 


t . 
1 


v t. 

-J-u.- 

u. i 

t. J+ 1 


(A-50) 


( u. , v. ) 

Vd+i Vi' 


A new line is considered which is tangent to the earth at 

and the process is repeated, resulting in a (j+l) (and better) 
estimate of aircraft location. When enough iterations have been taken, 
the aircraft's altitude and its location on the earth (subaircraft point) 
may be determined. 

The distance between the aircraft location estimate f u, , v„ 

1 A . A 


and the earth center (d g , 0) is Alfdg- u^ ) + v^ . Consequently 

'bll 

the aircraft's altitude at the j iteration is 


J 3 
r > 


h. = 
J 


1 C d 


V 2 1 

W) +V A. “ R e 
3 3 


(A-51) 


The subaircraft point /u^, , v^g^N is the intersection of the aircraft 

V 3 y 

location estimate-earth center line, given by 


v = 


d - u. 
e A. 
J 


u + 


d v. 
e A . 

i_ 

d - 

e A. 
0 


(A-52) 
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and the earth. 


u 


AE. 
J 


R u. + d h. 
e A e j 

J 

R + h. 
e J 


R V A 
e A , 


AE, 


R + h. 
e j 



(A- 53) 


This point is converted back to the nominal coordinate system, the 
x Q - y - z q system, by means of the coordinate transformation matrices, 

which are derived in Appendix B. 
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Appendix B 

COORDINATE TRANSFORMATION MATRICES 

1 - v y 0 - z 0 / x - y - 2 

The transformation from the nominal x Q - y 0 - z Q system to the 
body-fixed x-y-z system and vice versa are given by the following 
(see Fig. l) : 

I V 1 ?! \ 

V 2 ? 2 \ 

\ V 3 h ^3 





where 

= cos a 
= sin a sin y 
= sin a cos y 
= sin a sin g 

§2 = cos g cos y - cos a sin g sin y 

§2 = - cos g sin y - cos a sin g cos y 

71 ^ = - sin a cos g 

T) 2 = sin g cos y + cos a cos g sin y 

Tl^ = - sin g sin y + cos a cos g cos y 

oc, gj and y are the Euler angles involved. 


(B-2) 


(B-3) 
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2. x-y-z/u-v 

The purpose of this transformation is to facilitate the use of 
the iterative altitude determining procedure by performing all geometric 
manipulations in one plane rather than in three dimensions. Thus., this 
transformation results in rotating the x-y-z system in such a way that 
two axes are in the satellite-aircraft-earth center plane and one axis 
is perpendicular to it. 

The satellite-aircraft-earth center plane, P, is determined by 
those three points. Since the location of the aircraft is not known 
at this time, any point on the satellite-aircraft line will suffice. 

The equations for this line, derived in Appendix A, are 

x _ £ _ z 
a b c 

The point (a, b, c) is on the line. This point, together with the 
satellite, (0, 0, 0), and the earth center, (x e , y e , z e ), determine 
the equation for P in the following way: 



or (cy - bz ) x + (az - cx )y + (bx - ay )z = 0 (B-4) 

6 6 G G G G 

The intersection of the plane, P, with the x-y plane is given by 

x = m 1 y, (B-5) 

az - cx 

where m, = r— — • (B-6) 

1 bz - cy x ' 

e ^e 

Let 5 = tan ^m, (B-7) 

as shown in Fig. B-l. 

The z-axis is rotated an amount § such that the new y^axis is 
in P. The new x-, y-, and z-axes are called the x-, y-, and z-axes. 
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respectively. The matrix D provides for the proper rotation: 
/ * \ / * \ 




y 

z 


/ 


= D 


where 


0 \ 


0 


1 


where 


m. 


/ cos | -sin § 

D = sin § cos § 0 

A 0 

The intersection of P with the x - z plane is given by 
X = Mg Z, 


ay - bx 

e e 

l 2 _ (cy - bz ) cos 5 - (az - cx ) sin § 

G G 6 G 


(B-8) 


Let 


'll 


tan ^ nig, 


(B-9) 


(B-10) 


(B-ll) 


as shown in Fig. B-2. 


The y-axi£ is rotated an amount Tj such that the new z-axis is 
in P and the new x-axis is_perpendicular to P. The new x-, y-, and 
z-axes are called the x-, y- , and Z-axes, respectively. The matrix C 
provides for the proper rotation. 


/s\ 

y 

\ l ! 


= C 


/.\ 


y 

z 


where 
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c = 




cos T) 
0 

sin T| 


0 

1 

0 


-sin 7] 
0 

cos T] 


\ 


(B-12) 


The equation of the satellite-earth center line in the x-y-z 
coordinate system is given by 


where 


y = 


m„ = 


m^ z, 


x sin 5 + y cos 5 

0 G 


X cos l sin T] - y sin § sin Ij + z cos 7] 

0 0 .G 


(B-13) 


(B-14) 


Let 


v 


= tan 


-1 


m„ 


(B-15) 


as shown in Fig. B-3* 

The x-axis is rotated an amount v such that the new z-axis is 
along the satellite-earth center line, but with its positive direction 
pointed away from the earth. The new x-, y-, and z-axes are called the 

t t i , 

x- , y-, and z-axes, respectively. The proper rotation is provided for 
in the matrix, B. 


where 


B 


/ A 


1 =\ 

x' \ 


x \ 

y' 

= B 

y 

V'l 


l 1 / 


1 0 0 \ 

0 cos v -sin v ! 

y 0 sin v cos v j 


(B-l6) 


46 



Thus, the transformation from the x-y-z system to the x'-y'-z' 
system and vice versa is given by 


/ x '\ 


1 A 

y' 

= A 

y 

i z' 


i z 

\ i 


\ 


where 


and 


t*\ 


( x'\ 


T 


y 

= A 

r 

\ z 1 


\ z 'i 



r l 2 

r 13 ^ 

T 21 

i • 

r 22 

r 23 

\ r 31 

r 32 

r 33 1 


r ll 

=s 

COS 1 

cos 7] 



r 12 

= 

- sin § cos 

T| 


r 13 

= 

- sin 7] 



r 21 

= 

sin | 

cos V 

sin T] sin v 

r 22 

= 

cos l 

cosv + sin § 

sin T] sin v 

r 23 

= 

- cos 

7] sin 

V 


r 31 

= 

sin g 

sin v 

+ cos § 

sin T| cos v 

r 32 

= 

cos § 

sin v 

- sin c 

sin T| cos v 

r 33 

= 

cos T) 

cos V 




(B-17) 


(B-18) 


> (B-19) 
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The negative z'-axis is now called the q-axis (i.e., the u-axis 
coincides with the z'-axis but its positive direction is towards the 
earth), and the positive y'-axis is now called the v-axis, as shown 
in Fig. 4. In other words, 



(B-20) 


Consequently, the transformation from the u-v system to the 
x-y-z system is given by 


M 


M 


„T 


y 

= A 

V 

z 


l-u / 

\ ) 


\ / 


which reduces to 


r\ 


1 r 21 

r 31 \ 

y 

= ■ 

r 22 

i , . 

r 32 

\ z i 


\ r 23 

r 33 / 



(B-21) 


(B-22) 
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Appendix C 

TEST OF ITERATIVE ALTITUDE DETERMINING TECHNIQUE 


Consider Fig. C-l* which is the same as Fig. 4, except that 
it has additional angles marked. The values of At, u^g* and v^g will 
now be found* given 0* a ^* d g * R g and c (speed of light). 


d, cos(cp + a, ) = R sin a n 
1 Y 1' e 1 

d., sin(cp + a, ) + R cos a, = d 
,1 r e 1 e 

Combining these two equations yields 


(c-i) 

(C-2) 


, d cos a, - R 
, -1 e 1 e 

cp = tan — - t — — 

T d sin 


(C-3) 


From the law of sines* 
d 


sin(l8o°- 2cp) sin( 90 °- 0 + cp - a 1 ) 


(c-4) 


^ = ^l 

sin(0-9O°+ cp + c^) sin(90°- 0 + cp - a^) 

At can now be found using (C-3)* (C-4), and (C-5) from 

d l + *2 - d 
c 

2 (d g cos R g )(sin [a^+ 0]- cos cp) 
cos(a^+ 0 - cp) 


At = 


From the law of sines. 


_d 

sin a 


d 

e 

sin(l8o°- 0 - a) 


(C-5) 


(c-6) 


(c-7) 
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But from eq. (C-4) , 


sin a, sin a 1 


sin(l80°- 2cp) 
sin(90°- 0 + cp- 


(C-8) 


Again, from the law of sines. 


*2 


sin a. 


sin(90°+ 9) 


Combining (C-7)j (C-8), and (C-9) yields 


(C-9) 


a = tan 


-1 


2 sin 0 sin Qj sin cp 


cos(a i + 0 - cp) - 2 cos 0 sin 0^ sin cp 


(C-IO) 


From the figure. 


d e sin 0 = (R g + h) sin(0 + a) 


Therefore, 


d e sin 0 

h = sin (a + 0) ~ R e 


(c-11) 


(C-I2) 


Also, from the figure 


U AE = V R e cos 0 


(C-13) 


AE 


= R sin a 
e 


(C-14) 


Thus, considering the geometry of Fig. 4 (or Fig. C-l) and 
given 0, a-]_,d e , R e and c, it is possible to compute the corresponding 
values of At, h, u^ and v^. For instance, suppose 


0 = 8 . 0 ° 

0C ± = 58.5° 

d e = 42,237-92 km 
R g = 6371.26 km 
c = 2.997925 x 10 5 km/s 


(C-15) 


51 



Then, from equations (C-3), (C-6), (C-10), (C-12), (C-13), and (C-l4): 

At = 51-684414 (is 
h = 19-563931 km 
U AE = 38,946.876 km 

= 5455-454 km J 



Consider, now, the iterative altitude determining technique 
explained in Sec. 3-0* In this technique, the known quantities are 
At, d g , R g , and c. If these quantities have the same values as in 

eq. (c-15)above, and At has the value computed in (C-l6), the iterative 
technique yields 


h = 19.563931 km 
= 38,946.876 km ► 

V AE = 5455.454 km 


(C-17) 


It is noted that the values computed via the iterative technique in eq. 
(C-17) are the same as those computed in (C-l6). 


For a second example, let 


e = 6.22 0 



c 


39-5° 


7 same as above 



(C-l8) 


Using equations (C-3), (C-6), (C-10), (C-12), (C-13) and (C-l4) yields 
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(c-19) 


At 

h 


u 


AE 


AE 


48.522423 ms 
10.424577 km 
37,328.483 km 
4060.836 km 




With 0, d , R and c as above, and with At having the same value computed 
in (C-19), the iterative technique yields 


h = 10.424577 km 


AE 


AE 


37,328.483 km 
4060.836 km 


(C-20) 


Again, these values agree with those in (C-19). 

This test of the iterative technique has produced the same 
result in every one of the many situations considered. 
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